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FOREWORD 


Selective  calculations  have  been  made  for  burning  in  cracks  as  it  pertains 
to  solid  rocket  motor  propellants.  The  possibility  of  obtaining  pressures  high 
enough  to  cause  a  shock-to-detonation  transition  (SDT)  in  propellant  grains  is 
examined.  Variables  affecting  the  crack  combustion  process  which  were  selected 
for  study  are:  crack  shape,  location,  surface  roughness,  propellant  deformation, 
ignition  criterion,  and  burning  rate.  The  variables  are  evaluated  and  ordered 
in  groups  of  relative  importance.  Results  suggest  that  SDT  should  not  occur 
in  propellants  unless  the  granulation  of  the  grain  is  severe  enough  to  provide 
large  burning  surface  area. 
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INTRODUCTION 


Even  though  qualified  solid  rocket  motors  have  failed  due  to  fires,  case 
bursts  and  nozzle  expulsions,  no  service  motor  detonated  prior  to  19751.  The 
record  is  quite  impressive  considering  that  the  knowledge  of  solid  rocket 
motors  is  primarily  empirical.  One  reason  is  that  the  propellants  used  were 
of  moderate  performance  and  their  mechanical  strength  was  good. 

However,  propellants  under  development  since  1975  for  high  performance 
rocket  motors  contain  ingredients  which  are  more  detonable  including 
combustible  fillers  and  binders.  This  increases  the  general  sensitivity  of  the 
propellants  in  the  solid  content  (decrease  in  binder  content)  of  propellants 
which  can  result  in  poor  mechanical  properties  and  increased  sensitivity.  In 
fact,  several  motors  made  of  such  high  energy  propellants  detonated  high  order 
during  development. 

These  failures  have  been  carefully  investigated  and  there  are  several 
possible  scenarios  presented^: 

a.  Shear.  Damage  in  the  propel  1  ant- to-chamber  bond  region  creating  a 
bed  of  damaged  propellant  having  sufficient  granulation,  cleavage 

and  confinement  to  support  a  series  of  events  leading  to  deflagration- 
to-detonation  transition  (DDT). 

b.  Impact.  Detonation  of  expelled  propellant  fragments  on  nearby  objects 
leading  to  shock-to-detonation  transition  (SDT)  by  sympathetic  reaction 
from  these  expelled  fragment  reactions  to  the  propellant  bulk. 


1J.  F.  Kincaid  and  S.  J.  Jacobs,  "The  Fundamentals  of  High  Energy  Propellant 
Safety,"  Report  No.  10342.1  SP-20  CON.,  Dept,  of  Navy,  Strategic  Systems 
Project  Office,  1975. 

2 

J.  M.  Anderson,  J.  H.  Thacher  and  R.  B.  Wetherell ,  "Trident  Motor  Detonation 
Studies:  An  Overview,"  Presented  at  the  15th  JANNAF  Combustion  Meeting, 
Rhode  Island,  1978. 
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c.  Implosion. 

1.  Implosion-DDT,  wherein  the  grain  failed  near  the  case  bond  and 
is  deflected  inward  by  high  pressure  at  its  periphery,  causing 
breakup  or  cracking  in  the  centerport,  leading  to  DDT. 

2.  Implosion-SDT,  wherein  the  grain  collapse  is  sufficiently  rapid 
to  create  high  internal  shock  loads  at  the  centerport,  leading 
into  SDT. 

Given  these  scenarios,  attempts  have  been  made  to  analyze  them  theoretically 
and  in  laboratory  experiments  by  delineating  and  examining  the  essential  features 
in  detail.  Studies  include  DDT  in  granular  beds3,  DDT  and  SDT  in  both  cast 
and  pressed  explosives^,  cracks5,  and  impacts®. 

The  studies  also  include  the  investigation  of  convective  burning  processes 
in  propellant  cracks.  It  is  known  that  small  cracks  exposed  to  rocket  motor 
chamber  provide  more  surface  area  thereby  increasing  the  chamber  pressure. 

There  are  many  experimental  works  on  the  subject'. 

O 

Kuo  developed  a  one-dimensional  combustion  model  simulating  ignition 
and  combustion  in  solid  propellant  cracks  with  variable  cross  section  along 
the  length  of  the  crack.  Results  of  the  model  analysis  agreed  with  some 
experimental  evidence  at  low  pressures.  However,  although  these  studies 
concentrated  on  "flashdown",  a  flame  propagation  into  a  crack,  they  did  not 
consider  the  possibility  of  obtaining  high  enough  pressure  in  the  crack  to 
initiate  a  DDT  in  propellants. 

B.  D.  Hopkins,  N.  L.  Peterson,  M.  W.  Beckstead,  D.  T.  Pilcher  and  A.  G.  Butcher, 
"Mathematical  Modeling  of  the  DDT  Process,"  Presented  at  the  15th  JANNAF 
Combustion  Meeting,  Rhode  Island,  1978. 

A.  Forest,  "Numerical  Modeling  of  the  Deflagration-to-Detonation  Transition," 
Presented  at  the  15th  JANNAF  Combustion  Meeting,  Rhode  Island,  1978. 

®J.  L.  Prentice,  "Onset  of  Convective  Combustion:  Flashdown  in  Solid  Propellants," 
Presented  at  the  15th  JANNAF  Combustion  Meeting,  Rhode  Island,  1978. 

^R.  C.  Jensen,  "Investigation  of  High  Velocity  Impact  Initiated  Detonations," 
Presented  at  the  15th  JANNAF  Combustion  Meeting,  Rhode  Island,  1978. 

^A.  F.  Belyaev,  F.  K.  Bobolev,  A.  I.  Korotkov,  A.  A.  Sulimov,  and  G.  V.  Chuiko, 
Transition  From  Deflagration  to  Detonation  In  Condensed  Explosives,  Israel 
Program  for  Scientific  Translations,  Jerusalem,  1975. 

®K.  K.  Kuo,  A.  T.  Chen  and  T.  R.  Davis,  "Convective  Burning  in  Solid  Propellant 
Cracks,"  AIAAJ. ,  Vol .  16,  No.  6,  1978. 
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The  objective  of  this  report  is  to  evaluate  the  possibility  of  initiating 
high  order  detonation  in  propellants  through  crack  burning.  It  is  a  very 
complex  task  in  the  sense  that  no  set  of  definite  approaches  is  available, 
making  it  difficult  to  carry  out  systematic  experimental  study.  Furthermore, 
parameters  in  the  problem  (for  example,  crack  shapes  and  propellant  burning 
rate),  are  interdependent  and  add  more  complexity. 

In  this  report,  the  following  approaches  are  taken.  First,  factors  which 
affect  the  crack  burning  process  are  identified.  These  factors  should  include 
propellant  material  properties  as  well  as  crack  geometry.  Then  the  factors 
are  tested  for  their  relative  importance  in  the  pressurization  process  in 
burning  cracks.  For  this  purpose,  Kuo's  crack  combustion  model  was  modified  to 
include  the  effect  of  surface  roughness  and  propellant  deformation. 


IDENTIFICATION  OF  IMPORTANT  PARAMETERS 

Consider  a  simplified  crack  combustion  system  shown  in  Fig.  1.  An  amplified 
view  of  a  crack  is  shown  to  face  the  rocket  chamber  filled  with  hot,  high 
pressure  gas.  The  hot  gas  initially  penetrates  into  the  crack  cavity,  heating 
the  propellant  surface  through  convective  heat  transfer.  Some  of  the  heat  is 
lost  through  conductive  heat  transfer  into  the  propellant  body.  If  the  convective 
heat  transfer  is  higher  than  the  heat  loss  through  conductive  heat  transfer, 
the  surface  temperature  increases  until  it  reaches  the  ignition  temperature  of 
the  propellant.  When  ignited,  the  propellant  burns  at  a  rate  dependent  on  the 
local  gas  pressure  whose  rate  of  buildup  depends  on  the  burning  rate  and  the 
geometry  of  the  crack.  A  rough  surface  provides  a  large  burning  area  and  hence 
a  large  gasification  rate.  As  the  burning  rate  increases,  a  compressive  wave 
is  formed  which  enters  and  propagates  along  the  crack  ahead  of  the  flame  front 
(ignition  front).  The  crack  cavity  is  further  pressurized  by  the  gasification 
of  the  solid  propellant  and  by  reflection  of  the  compression  wave  at  the  end 
of  the  crack.  The  pressurization  continues  until  the  crack  width  is  sufficiently 
enlarged  by  consumption  of  propellant  so  that  the  gas  escape  rate  through  the 
crack  opening  is  larger  than  the  gas  generation  rate  inside  the  crack.  In  the 
current  discussion,  crack  propagation  (increase  of  crack  length)  is  not 
considered. 

The  crack  combustion  system  can  be  analyzed  by  solving  a  set  of  partial 
differential  equations  (conservation  equations  for  gas  phase)  with  appropriate 
initial  and  boundary  conditions,  containing  non-homogeneous  mass,  momentum  and 
energy  source  terms.  Parameters  which  influence  the  crack  combustion  process 
are: 


a.  Geometry.  There  are  several  independent  parameters  describing  the 
geometry  of  the  crack.  These  are  the  crack  length,  L,  the  crack 
width,  W,  and  the  roughness  of  the  crack's  inner  surface.  Small 
cracks  may  branch  off  the  main  crack  to  provide  additional  surface 
area  for  combustion.  Also,  the  deformation  of  the  propellant  due  to 
the  gas  pressure  in  the  crack  plays  an  important  role  in  the  crack 
combustion  process. 
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b.  Initial  Conditions  (Tn,p0)  In  most  cases,  initial  conditions  inside 
the  crack  are  not  important  to  the  crack  combustion  process.  The 
initial  pressure  is  small  compared  to  the  chamber  pressure  and  does 
not  present  significant  resistance  to  the  hot  chamber  gas  penetration 
into  the  crack. 

c.  Boundary  Conditions.  The  boundary  conditions  are: 

1.  The  pressure  at  the  inlet  of  the  crack  equals  the  chamber  pressure 
or  there  exists  a  choked  flow  condition  at  the  crack  inlet. 

2.  The  gas  velocity  at  the  closed  crack's  end  is  zero. 


CHAMBER 

> 

INITIAL  HOT  GAS 
PENETRATION 


FIGURE  1  SCHEMATIC  CRACK  BURNING  PROCESS 


If  the  crack  is  initially  exposed  to  the  rocket  chamber  when  the  rocket 
motor  is  being  ignited,  then  the  chamber  pressure  (crack  inlet  boundary 
condition)  is  a  slowly  increasing  function  of  time  until  it  reaches  a  motor 
operating  pressure  (arbitrarily  chosen  here  to  be  about  10  MPa).  The 
characteristic  time  scale  is  of  the  order  of  0.1  seconds.  Kuo°  has  shown  that 
under  this  condition,  the  pressure  in  the  crack  is  almost  identical  to  the 
chamber  pressure,  indicating  that  there  is  no  additional  pressure  buildup  due 
to  the  presence  of  the  crack.  For  this  reason,  the  case  where  the  crack  is 
initially  exposed  to  the  rocket  chamber  is  not  discussed  here.  The  cracks 
considered  are  those  suddenly  opened  to  the  already  established  chamber  pressure 
of  10  MPa  (embedded  cracks  with  one  side  of  the  crack  being  exposed  to  rocket 
chamber  pressure  by  burning  of  the  propellant)  or  -  suddenly  created  due  to 
mechanical  failure.  The  boundary  condition  (1)  stated  above  is  modified  in 
this  casd  as  follows: 

1.  The  pressure  at  the  inlet  of  the  crack  is  10  MPa  or  there  exists 
a  choked  flow  condition  at  the  crack  inlet. 


8 


NSWC  TR  79-485 


d.  Ignition  and  Burning  Characteristics  of  Propellants.  Burning 
characteristics  of  the  propellants  also  contributes  to  the 
pressurization  process  as  in  solid  propellant  rocket  motors  in 
which  burning  rate  and  rocket  motor  geometry  together  establish 
a  stable  operating  chamber  pressure.  Although  ignition 
characteristics  are  not  important  in  rocket  motors  because  the 
ignition  delay  time  is  relatively  short  compared  to  the  burning 
time  to  maximum  pressure,  they  will  be  considered  in  this  paper 
because  the  characteristic  time  in  the  crack  combustion  process 
may  be  of  the  same  order  of  magnitude  as  the  ignition  delay  time. 

In  the  following  treatment,  only  items  (a)  and  (d)  of  the  above  will  be 
studied.  As  discussed  in  the  INTRODUCTION,  it  is  impossible  to  evaluate  all 
parameters  in  a  comprehensive  way  because  of  the  nonliniarities  of  the  problem 
and  the  interdependency  of  the  parameters.  The  approach  chosen  in  this  paper 
is  to  examine  each  parameter  singly  and  establish  its  importance  with  a  modifi¬ 
cation  of  Kuo's  Crack  Combustion  Code.  The  code  is  explained  in  detail  in  the 
Appendix. 


EVALUATION  OF  PARAMETERS 


BASELINE  CASE 


Initial  crack  shape  was  studied  first  disregarding  propellant  deformation. 
Two  cases  are  considered  for  the  geometrical  effects.  Parameters  used  are 
listed  in  Table  1.  The  length  to  width  ratio  of  the  crack,  L/W,  is  about 
300  in  Case  I,  and  1000  in  Case  II. 

Case  I.  The  details  of  the  crack  combustion  process  computed  by  Kuo's 
model  are  illustrated  in  Figs.  2-5  for  Case  I  and  compared  with  Fig.  6  for 
Case  II.  Fig.  2  shows  a  p-x  diagram  at  various  times  for  Case  I.  x  =  0  at 
the  crack  entrance  and  x  =  L  =  15.05  cm  at  the  crack  tip.  As  stated 
previously,  the  chamber  pressure  is  maintained  at  100  bars  (10  MPa).  Note 
that  the  initial  pressure  wave  travels  for  about  150  usee  before  it  is  reflected 
at  the  crack  tip.  Dashes  on  the  p-x  profiles  for  50  ~  150  usee  indicate  the 
position  of  flame  front  (or  ignition  front)  at  these  times.  The  figure  therefore 
shows  that  the  ignition  front  follows  the  pressure  wave  front  closely  and,  at 
about  150  usee,  the  entire  crack  is  ignited.  The  pressure  at  the  time  is 
more  or  less  constant  throughout  the  length  of  the  crack  but  steadily  increases 
as  the  wave  is  reflected  from  the  crack  tip.  Even 

after  a  rarefaction  wave  comes  into  the  crack  after  300  usee,  this  pressurization 
continues.  Pressurization  will  continue  until  the  crack  width  is  sufficiently 
enlarged  for  the  mass  generation  rate  in  the  crack  to  equal  the  rate  of  mass 
escape  out  of  the  crack.  The  pressure  peak  will  be  discussed  further  when 
propellant  deformation  is  studied. 
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TABLE  1  INPUT  PARAMETERS  USED  FOR  SIMULATION 
W  =  crack  width  =  0.0508  cm  for  Case  I,  and  =  0.015  cm  for  Case  II 
L  =  crack  length  =  15.05  cm 
Bx  =  Body  force  =0.0  g^/g 

3 

b  =  co-volume  =1.0  cm  /g 

Mw  =  molecular  weight  =  23.97  g/gr-mole 

y  =  specific  heat  ratio  =  1.233 

3 

Ppr  =  propellant  density  =  1.741  g/cm 

Apr  =  thermal  conductivity  =  0.55  x  10"3  cal/sec-cm-°K 

p 

a  =  thermal  diffusitivity  =  0.001  cm  /sec 

r  * 

er  =  relative  roughness 

=  surface  roughness/hydraulic  diameter  of  the  crack  =  0.001 
Tf  =  flame  temperature  =  3000°K 
Tcri  =  critical  temperature  =  780°K 
^ign  =  i9n''tion  temperature  =  1000°K 
T|  =  initial  temperature  =  298°K 

a  =  pre-exponential  factor  .  2  n 

in  burning  rate  law  =  2.752  x  10"^  (c.m/sec)/(g/cm  ) 

n  =  exponent  =  0.683 

Kg  =  erosive  burning  constant  =  15.9  cm3  -  °K/cal 
B  =  erosive  burning  exponent  =  125 
Tc  =  chamber  temperature  =  2800°K 
Pc  =  chamber  pressure  =  100  bars  for  t  >  0 


10 


100  bars 


0  5  10  15 


x,  cm 

FIGURE  2  CALCULATED  PRESSURE  PROFILES 
AT  VARIOUS  TIME  —  BASELINE 
CASE  I.  IGNITION  FRONT  IS  SHOWN 
AT  BOTTOM  OF  THE  FIGURE. 
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400  /k  sec 
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Figure  3  shows  the  gas  temperature  profile.  It  increases  from  the 
initial  value  of  298°K  to  around  3000°K  as  the  first  compression  wave 
propagates  into  the  crack  and  the  temperature  remains  approximately  constant 
throughout  the  process.  One  exception  of  this  is  that  at  the  end  of  the 
crack,  the  cold  gas  which  was  initially  in  the  crack  is  less  than  3000°K 
although  it  is  compressed  at  the  end  of  the  crack  to  form  a  cool  pocket  in 
that  region. 


FIGURE  3  TEMPERATURE  PROFILES  AT  VARIOUS 
TIMES  —  BASELINE  CASE  I 


Gas  velocity  profiles  as  a  function  of  space  at  various  times  are  depicted 
in  Fig.  4.  The  initial  compression  wave  traveling  inward  is  evident  from  the 
positive  velocity  profiles  at  50,  100  and  150  psec.  At  150  psec,  some  part 
of  the  gas  near  the  open  end  escapes  to  the  chamber.  At  300  psec,  the  velocity 
of  the  bulk  of  the  gas  is  almost  at  rest  (in  this  particular  case,  slightly 
positive)  and  near  the  open  end  the  gas  is  rushing  out.  From  then  on,  the 
bulk  of  the  gas  maintains  accelerating  gas  profiles  toward  the  open  end. 

The  change  of  the  crack  width  in  time  is  shown  in  Fig.  5.  The  original 
crack  width  is  0.0508  cm.  The  width  increases  fastest  near  the  open  end 
because  that  end  is  ignited  first  and  because  high  gas  velocity  increases 
the  erosive  burning  effect  there  (Refer  to  Fig.  4). 
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FIGURE  4  VELOCITY  PROFILES  —  BASELINE  CASE  I 
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FIGURE  5  CRACK  WIDTH  —  BASELINE  CASE  I 


Case  II.  Case  II  results  are  shown  in  Fig.  6  for  comparison.  L/W  ratio 
for  Case  II  is  1000  as  compared  to  300  for  Case  I.  The  figure  displays  P-x 
profiles  in  the  crack  at  various  times  as  did  Fig.  2.  The  smaller  crack 
width  (or  the  larger  L/W)  -  increases  the  pressure  in  the  crack  significantly 
by  the  time  the  ignition  first  reached  the  crack's  tip,  the  peak  pressure  has 
already  exceeded  1  kbar  and  is  rising  rapidly  (at  225  usee,  the  peak  pressure 
is  around  8  kbar).  This  illustrates  the  importance  of  the  crack  geometry  in 
the  pressurization  of  the  crack. 
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IGNITION  CRITERION 

Results  shown  thus  far  indicate  that  the  ignition  front  (or  the  flame 
front)  closely  follows  the  first  compression  wave  propagating  into  the  crack. 
Since  the  major  pressurization  in  cracks  occurs  after  the  entire  crack  is 
ignited,  moderate  changes  in  ignition  phase,  or  in  ignition  characteristics 
are  not  expected  to  change  the  crack  combustion  process  significantly.  This 
is  shown  in  Fig.  7.  The  situation  is  identical  to  the  baseline  Case  I  except 
that  the  ignition  temperature  is  lowered  by  100°K.  The  trend  is  similar  to 
the  baseline  case  except  that  due  to  the  choice  of  a  lower  ignition  temperature, 
the  first  compression  front  propagates  a  little  faster  than  in  the  baseline 
case.  There  is  roughly  a  25  ysec  gain  which  persists  throughout  the  combustion 
process. 


SIMILARITY 

There  has  been  an  emphasis  on  the  ratio  L/W  rather  than  on  L  or  W 
separately  in  crack  combustion  studies9.  When  the  pressurization  in  the 
crack  is  moderate  at  low  L/W,  the  propellant  deformation  is  not  severe,  and 
as  the  first  compression  wave  into  the  crack  is  reflected  from  the  crack  tip, 
the  crack  is  nearly  uniformly  pressurized  and  ignition  is  completed,  the 
following  processes  depend  entirely  on  fluid  mechanics  with  the  mass  generation 
term  dependent  upon  fluid  dynamic  variables  only.  The  resulting  partial 
differential  equation  describing  the  flow  process  from  that  point  on  can  be 
scaled  with  L/W  as  long  as  thermodynamic  variables  in  the  crack  are  approximately 
constant  when  the  ignition  is  completed.  This  deduction  has  been  verified  by 
examining  Fig.  8  which  is  identical  with  Case  I  except  that  L  and  W  are  both 
3/10  times  smaller. 

However,  at  high  L/W,  the  uniform  initial  condition  cannot  be  assumed 
after  the  compression  wave  is  reflected  at  the  crack  tip  as  one  observes 
a  severe  pressure  gradient  in  the  crack.  Furthermore,  propellant  deformation 
plays  a  more  important  role  in  this  case  for  a  similarity  to  hold.  Therefore, 
at  high  L/W  ratio,  L  and  W  should  be  independently  studied  rather  than  a 
single  parameter  representation  of  L/W. 

PROPELLANT  DEFORMATION 

The  crack  geometry  is  dependent  upon  many  factors.  The  initial  geometry 
is  important,  and  as  burning  progresses,  the  increased  width  due  to  propellant 
consumption  also  plays  a  role.  Another  important  factor  is  the  increased 
width  caused  by  high  internal  pressure  deforming  the  propellant  body 
surrounding  the  crack.  The  amount  of  the  propellant  deformation,  or  the 
enlargement  of  the  crack,  for  a  given  propellant  depends  on  the  surrounding 
geometry.  Consider  a  geometry  shown  in  Fig.  9. 

g 

A.  T.  Belyaev,  A.  I.  Korotkov,  A.  A.  Kulimov,  M.  K.  Sukoyan  and  A.  V.  Obmenin, 
"Development  of  Combustion  in  an  Isolated  Pore,"  Combustion,  Explosion  and 
Shock  Waves ,  Vol .  5. 
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FIGURE  7  PRESSURE  PROFILES  —  LOWER  IGNITION  TEMPERATURE 
THAN  THE  BASELINE  CASE  I 
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FIGURE  8  PRESSURE  PROFILES  —  SIMILAR-SHAPED  CRACK 

AS  THE  BASELINE  CASE  I  BUT  3/10  TIMES  SMALLER 
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A  crack  of  length,  L,  and  width,  W,  is  sandwiched  between  two  propellant 
slabs  of  thickness  of  Z/2  each  which  are  backed  by  rigid  walls.  If  Z  is 
small  enough,  for  example,  of  the  order  of  1  cm,  so  that  the  wave  in  the 
propellant  propagates  laterally  in  a  much  shorter  time  than  the  characteristic 
time  of  the  crack  combustion,  a  quasi-equilibrium  approximation  can  be 
introduced  in  calculating  propellant  deformation.  With  a  propellant  bulk 
modulus,  K,  neglecting  any  viscoelastic  effect,  one  can  estimate  the  change 
in  the  propellant  thickness, 

AZ  =  -  l  P 
IT 

where  p  is  the  pressure  in  the  crack.  The  Kuo  code  was  modified  to  include 
this  effect. 

Normally,  in  large  rocket  motors  it  is  difficult  to  achieve  a  condition 
such  as  shown  in  Fig.  9,  with  a  single  crack  and  l  will  be  normally  larger 
than  1  cm.  As  l  increases,  AZ  becomes  larger  thus  making  the  width  of  the 
crack  larger.  This  reduces  the  magnitude  of  the  peak  pressure  achievable  in 
the  crack  combustion.  However,  the  above  condition  can  be  qualitatively 
achieved  in  parallel  cracks  of  equal  lengths  and  equal  widths. 


FIGURE  9  IDEALIZED  SINGLE  CRACK  GEOMETRY  IN  MULTIPLE 
CRACK  SITUATION 

Since  a  large  l  decreases  the  maximum  pressure  in  the  crack,  for  the 
purpose  of  studying  the  possibility  of  obtaining  high  pressure  for  shock 
initiation  of  detonation  in  the  propellants,  the  smallest  value  of  l  was 
examined.  With  the  crack  length,  L  =  15.05  cm  and  the  width,  W  =  150  pm,  a 
value  of  0.5  cm  for  l  was  investigated.  With  a  solid  particle  size  of  the 
order  of  200  pm,  it  is  difficult  to  conceive  that  any  propellant  thickness, 

1%  less  than  0.5  cm  between  many  parallel  cracks  can  be  structurally  strong. 

Weak  propellants  will  granulate.  Once  the  propellant  is  granulated,  crack 
combustion  model  is  no  longer  applicable,  and  convective  combustion  in  granular, 
porous  bed  should  be  studied  instead.  Convective  combustion  in  a  granular  bed 
is  beyond  the  scope  of  this  report. 
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The  result  for  the  case,  after  modification  of  Kuo's  model,  l  =  0.5  cm, 

K  =  600  Kpsi  is  shown  in  Fig.  10.  It  can  be  compared  with  Fig.  6  (no  deforma¬ 
tion).  Propagation  of  the  initial  compression  wave  into  the  crack  is  almost 
identical  in  both  cases  at  early  times  due  to  the  relatively  low  pressure  (small 
propellant  deformation).  The  ignition  of  the  propellant  is  achieved  at  the 
leading  edge  of  the  compression  wave  as  before. 

At  later  times  (>175  ysec),  the  cases  differ.  In  the  non-deforming  case 
(Fig.  6),  the  pressurization  is  extremely  rapid  and  a  slow-down  of  the  pressuri¬ 
zation  is  not  observed.  In  the  deforming  case  (Fig.  10),  a  maximum  pressure 
is  observed  near  400  ysec,  after  which  the  pressure  decreases.  The  peak 
pressure  (at  400  ysec)  is  about  5  kbars.  If  this  pressure  is  maintained  long 
enough,  some  propellants  may  react  chemically  at  hot  spots. 

ROUGH  SURFACE,  HIGH  BURNING  RATES 

Thus  far,  the  study  has  concentrated  on  cracks  of  relatively  smooth  surface. 
However,  in  view  of  the  small  crack  width  (of  the  order  of  magnitude  of  some 
of  the  larger  particles  in  the  propellant),  it  is  expected  that  any  given 
crack  would  have  rough  surfaces.  The  effect  of  the  rough  surface  is  twofold: 
it  retards  the  flow  through  large  drag  forces  and  it  forms  a  larger  burning 
surface  area  once  the  propellant  is  ignited.  Many  small  cracks  branching 
sideways  from  the  main  crack  may  also  exist.  These  would  certainly  increase 
the  available  burning  surface  area  and  thereby  increase  pressurization  in  the 
main  crack. 

Exact  analysis  of  the  above  situation  would  be  difficult  and  time  consuming. 
Instead,  a  simplistic  approach  was  taken  to  account  for  the  effects  of  surface 
roughness  and  many  small  branching  cracks.  To  simulate  the  situation,  the 
burning  rate  was  increased  by  an  arbitrary  factor,  a,  while  the  regression 
rate  was  kept  as  before.  This  would  have  the  same  effect  as  having  a  times 
larger  burning  surface  area  with  the  same  crack  width.  Therefore,  the  result 
of  this  computation  would  be  expected  to  produce  more  severe  pressurization 
than  that  of  just  higher  burning  rate  with  correspondingly  higher  regression 
rate  (wider  crack  would  result  from  faster  burning).  Since  here  the  major 
concern  is  on  the  worst  possible  condition  achievable  in  propellant  crack 
burning,  the  latter  can  be  ignored  in  this  study. 

The  results  of  computations  for  values  of  a  =  2  and  a  =  3  are  shown  in 
Fig.  11  and  Fig.  12.  Compared  with  Fig.  10  (a  =  1),  it  may  be  noted  that  the 
maximum  pressure  has  a  roughly  linear  functional  relationship  to  values  of  a, 
for  the  range  considered.  A  simple  mass  balance  equation  (gas  generation  due 
to  burning  vs  gas  escape  through  the  crack  opening)  also  confirms  this 
rough  linear  relationship  (not  shown  here). 

The  pressure  peak  observed  in  Fig.  12  is  high  enough  and  has  a  long  enough 
duration  to  cause  some  propellants  to  transit  to  detonation.  A  higher  a 
value  would  worsen  the  situation.  However,  this  higher  value  of  a  means 
very  extensive  cross-cracking  in  a  multiple  crack  area.  The  analysis  performed 
here  would  become  less  effective  for  this  case  and  probably  a  convective  burning 
analysis  in  granular  beds  would  be  better  employed  for  this  purpose. 


20 


NSWC  TR  79-485 


CONCLUSIONS 

Fundamental  mechanisms  in  crack  combustion  processes  were  examined  to 
determine  the  possibility  of  obtaining  high  enough  pressure  to  cause  a 
transition  to  detonation.  Several  parameters  were  tested  for  their  effects  on 
the  pressurization  in  cracks.  They  are  summarized  as  follows. 

1.  The  crack  has  to  be  an  embedded  one  or  a  suddenly  created  one  to 
generate  high  pressure  in  the  crack. 

2.  The  ratio  of  the  crack  length  to  width  is  an  important  parameter. 

At  higher  values  of  this  ratio  the  length  and  the  width  must  be 
considered  separately. 

3.  The  propellant  thickness,  lt  between  multiple  parallel  cracks  is 
important. 

4.  The  factor  a,  which  reflects  the  amount  of  surface  area  available 
for  combustion  in  cracks  shows  an  almost  linear  relationship  with 
the  peak  pressure  attainable  in  the  crack  in  a  multiple  parallel 
network. 

5.  Among  propellants  with  the  same  exponent  for  burning  rate  expression, 
propellants  with  larger  branching  cracks  with  lower  absolute  burning 
rate  would  generate  higher  pressure  in  the  cracks  than  propellants 
of  smooth  surface  with  higher  burning  rate. 

6.  Ignition  criterion  is  the  least  important  of  the  parameters  studied. 

7.  Among  the  parameters  discussed  above,  some  are  limited  in  their  values 
to  be  considered  realistic  in  the  current  study.  The  propellant 
thickness,  l,  cannot  become  too  small,  nor  can  the  degree  of  the 
crack  be  too  extensive;  otherwise  the  situation  no  longer  represents 

a  single  crack  combustion  process.  Rather,  it  must  be  studied  as  a 
granular  bed  combustion  process. 

8.  Peak  pressure  high  enough  to  cause  propellants  to  transit  to  detonation 
has  not  been  observed  in  isolated  smooth,  short  cracks,  but  only  in 
long,  rough  cracks  in  a  multiple  crack  situation  with  very  thin 
propellant  layers  in  between  cracks. 

9.  Finally,  experimental  verification  of  the  predictions  is  desired  in 
order  to  test  the  assumptions  used  in  this  model  and  to  provide  more 
definitive  directions  for  future  studies.  Also,  information  on  the 
nature  of  real  motor  cracks  should  be  provided  to  guide  further 
investigation. 
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APPENDIX  A 

CRACK  COMBUSTION  CODE  (KUO) 


Kuo's  crack  combustion  code  is  described  in  a  simplified  form  in  the  following. 
Basic  Assumptions 

The  following  basic  assumptions  are  used  in  the  theoretical  model: 

1.  All  chemical  reactions  occur  near  the  propellant  crack  surface,  in  a 
combustion  zone  that  is  so  thin  that  it  can  be  considered  as  a  plane.  That  is, 
the  thickness  of  the  combustion  zone  above  the  propellant  surface  is  much 
smaller  than  the  gap  width  in  the  propellant  crack. 

2.  Rate  processes  at  the  propellant  surface  are  quasisteady  in  the  sense 
that  the  characteristic  times  associated  with  the  gaseous  flame  and  preheated 
propellant  are  short  compared  to  that  of  the  pressure  transient  variation. 

3.  The  gases  flowing  in  the  propellant  crack  obey  the  Calusius  or  Noble-Abel 
equation  of  state.  This  dense  gas  relation  can  adequately  describe  the  departure 
from  the  ideal  gas  law  at  high  pressures. 

4.  The  bulk  flow  in  the  pore  is  one-dimensional.  Recent  experiments 
conducted  by  Kuo  et  al.  *  have  revealed  that  the  one-dimensional  assumption 
is  appropriate  for  flow  in  cracks  when  the  pressure  gradient  along  the  crack 
is  greater  than  a  critical  value,  which  normally  is  extremely  small. 

5.  The  compressibility  of  the  solid  propellant  is  neglected.  The  change 
of  port  area  in  propellant  cracks  is  solely  due  to  the  propellant  surface 
regression. 

Theoretical  Modeling 

To  describe  the  gas-phase  behavior  inside  a  solid-propellant  crack,  the  mass, 
momentum,  and  energy  equation  in  unsteady,  quasi-one-dimensional  forms  have 
been  developed  based  upon  the  balance  of  fluxes  in  a  control  volume  within 
the  propellant  crack. 

The  mass  conservation  equation  is 

3 (pA  )  3  (puA  ) 

3t  +  3x  =  rbpprpb  ^ 


A-l 
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The  momentum  conservation  equation  is 


L_  (puA  )  +  —  (PAnu2)  =  -  A„  9JP 
3t  P 


-  p.. T  COS0, , 

rw  w  w 


P  Pax 

(pprrbPb)Vgfsin6w 


(2) 


The  energy  conservation  equation  written  in  terms  of  the  total  stored 
energy  (internal  and  kinetic)  per  unit  mass,  E,  is 

ft  <PV>  +fl  <PVE)  " 

-  h  (VU)  +  pprVbhf 

-Vb(T-Tps>  (3) 

After  rearranging  the  conservation  equations,  a  set  of  three  governing 
equations  is  obtained,  i.e.,  velocity-variation,  temperature-variation,  and 
pressure-variation  equations.  These  three  governing  equations  are  first-order, 
coupled,  inhomogeneous,  and  nonlinear  partial  differential  equations.  The 
unknowns  that  vary  strongly  with  respect  to  the  spatial  variables  are  gas 
temperature,  pressure,  velocity,  and  the  temperatures  at  the  propellant 
surface  and  inert  wall.  There  are  other  unknown  parameters  that  vary  weakly 
with  respect  to  spatial  variables.  These  parameters  include  the  port  area  of 
the  crack,  the  burning  rate  of  the  propellant,  the  convective  heat  transfer 
and  friction  coefficients,  the  gas  flow  angle,  the  local  blowing  velocity 
normal  to  the  propellant  surface,  and  the  burning  perimeter  at  various 
locations  along  the  crack. 

To  solve  all  the  unknowns  mentioned  above,  additional  equations  and 
boundary  conditions  together  with  some  physical  input  functions  must  be  specified. 
The  propellant  surface  temperature  at  a  fixed  location  along  the  crack  before 
the  attainment  of  ignition  is  calculated  from  the  solid-phase  heat  conduction 
equation  which  is  written  in  the  unsteady  one-dimensional  form 

3Tpr(t.y)  32Tor(t.y)  (4) 

at  "  V  g y2 

where  the  length  variable  y  is  measured  perpendicular  to  the  local  propellant 
crack  surface.  The  initial  and  boundary  conditions  are 

Tpr(°.y)  -  Tpl.  (5) 


A- 2 
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(6) 


3Tor  hc(t)  , 

■Sf  (t-°>  *-t£-  fT<‘>  -  TPs(t)}  <7> 

The  heat  conduction  equation  is  solved  either  by  using  an  integral  method 
employing  a  third-order  polynomial  or  by  direct  numerical  solution  of  Eqs.  (4-7) 
with  variable  mesh  size  in  the  subsurface. 

For  the  gas  phase,  the  Nobel-Abel  equation  is  used  for  the  equation  of 
state: 


P ( 1/p-b )  =  RT  (8) 

For  the  solid  phase,  the  statement  of  a  constant  density  for  the  solid  propellant 
serves  as  the  equation  of  state. 

The  expression  for  the  local  convective  heat- transfer  coefficient,  hc,  is 
deduced  from  the  conventional  Dittus-Boelter  correlation  for  turbulent  flow 
in  pipes. A_z  Variation  of  the  physical  properties  of  the  gas  across  the  boundary 
layer  is  considered  by  evaluating  the  properties  at  an  average  film  temperature, 
Tflf.  Before_local_ignition  hcp  =  hew  =  hc,  whereas  after  local  ignition 
Hep  =  0  and  hCw  =  hc.  The  Prandtl  number  used  in  the  Dittus-Boelter  correlation 
is  calculated  from  Svehla's  equation. The  correlation  for  the  friction 
coefficient  used  in  this  study  is  the  modified  form  of  the  well-known  Coel brook 
formula,  4  which  is  discussed  in  Reference  A-2.  Zero  wall  friction  is  assumed 
at  the  burning  surface  inside  the  crack  due  to  large  friction  attenuation 
caused  by  the  surface  blowing. 

In  the  absence  of  a  more  adequate  erosive  burning-rate  expression,  the 
Lenoir-RobillardA_5  burning-rate  formula  is  used  to  account  for  erosive  burning 
effect.  After  Tps  is  solved,  a  simplified  ignition  temperature  criterion  is 
used  to  determine  the  burning  condition  of  the  solid  propellant  along  the  crack. 
To  avoid  the  step  function  change  in  burning  rate,  a  two-temperature  criterion 
is  used  to  achieve  full  ignition  with  a  finite  time  interval;  when  Jps  Is,  equal 
to  TCy.-j,  the  propellant  starts  to  gasify.  As  soon  as  Tps  reaches  T-,- gn ,  the 
full  ignition  condition  is  reached. 

The  initial  conditions  required  for  the  solution  of  the  system  of  governing 
equations  are 


u(0,x)  =  uj  T(0,x)  =  li  P(0,x)  =  Pi  (9) 

where  Ui ,  Tj ,  and  Pj  are  the  initial  velocity,  temperature,  and  pressure  in 
the  propellant  crack. 
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The  number  of  boundary  conditions  required  depends  upon  the  flow  direction 
and  Mach  number  at  the  opening  of  the  crack.  When  the  gas  in  the  rocket  chamber 
flows  subsonically  into  the  propellant  crack,  two  boundary  conditions  at  the 
opening  of  the  propellant  crack  are  specified  to  represent  the  interface 
conditions  at  the  crack  entrance;  and  one  boundary  condition  at  the  closed 
end  of  the  propellant  crack  is  specified  to  indicate  that  the  gas  velocity 
is  zero.  In  the  mathematical  form,  they  are  represented  as  follows: 

P(t,0)  =  Pc(t)  T(t,0)  =  Tc(t)  u(t,xL)  =  0  (10) 

When  the  flow  of  the  gas  out  of  the  crack  is  subsonic,  two  boundary 
conditions  are  specified.  They  are  written  as 

P(t,0)  =  Pc(t)  u(t,xL)  =  0  (11) 

When  the  outflowing  gas  is  supersonic,  only  one  boundary  condition  can  be 
specified,  that  is 


u(t,xL)  =  0  (12) 

To  solve  the  system  of  governing  equations  with  a  second-order  numerical 
scheme,  extraneous  boundary  conditions  are  required  in  addition  to  the  physical 
ones.  The  procedure  to  attain  these  conditions  is  discussed  in  the  following 
section 

Numerical  Solution  Technique 

The  governing  equations,  describing  the  flame-spreading  and  combustion 
processes  in  a  propellant  crack,  were  found  to  be  totally  hyperbol ic.A~° 

The  numerical  techniques  developed  were  chosen  on  the  basis  of  stability,  maximum 
accuracy,  and  computational  efficiency.  Recently  obtained  experience  in  solving 
hyperbolid  partial  differential  equations^"2  was  utilized,  and  a  generalized 
implicit  scheme^-'  based  on  central  differences  in  spacewise  derivatives  was 
chosen  to  solve  the  governing  equations.  A  quasi  linearization  method,  used  to 
linearize  the  inhomogeneous  terms  in  the  governing  equation,  was  combined  with 
a  stable  predictor-corrector  procedure. A_2  Finally,  the  governing  equations,  in 
their  finite-difference  form,  are  arranged  in  a  block-tridiagonal  matrix  form,A~2 
which  allows  an  efficient  computation. 

The  use  of  central  difference  formulation  for  all  spacewise  derivatives 
in  the  governing  equations  require  six  boundary  conditions  for  the  solution. 

This  implies  that,  when  the  gas  flows  into  the  crack,  three  extraneous  boundary 
conditions  are  required  in  addition  to  the  three  physical  ones  described 
by  Eq.  (10).  The  three  extraneous  boundary  conditions  are  used  to  determine 
the  flow  velocity  at  the  opening  of  the  crack,  and  the  pressure  and  temperature 
at  the  closed  end  of  the  crack. 
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When  the  gas  flows  out  of  the  crack  subsonically,  in  addition  to  the  two 
physical  boundary  conditions  described  in  Eq.  (11),  two  extraneous  boundary 
conditions  are  required  at  the  opening  of  the  crack  for  the  determination  of 
the  pressure  and  temperature.  Also,  two  extraneous  boundary  conditions  are 
required  at  the  closed  end  of  the  crack. 

When  the  outflowing  gas  is  supersonic,  in  addition  to  the  physical  boundary 
condition  given  by  Eq.  (12),  three  extraneous  boundary  conditions  are  needed 
at  the  opening  of  the  crack  and  two  other  extraneous  boundary  conditions  are 
needed  at  the  tip  of  the  crack. 

These  extraneous  boundary  conditions  for  the  hyperbolic  system  of  governing 
equations  are  derived  from  the  compatibility  relations  at  the  boundaries.  The 
relations  are  obtained  by  transforming  the  governing  equations  using  the  method 
of  characteristics. ^"8  jhe  compatibility  relations  together  with  the  physical 
boundary  condition  form  a  closed  system  to  determine  the  velocity,  temperature, 
and  pressure  at  both  ends  of  the  crack. 
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NOMENCLATURE 

p 

Ap  =  cross-sectional  area  of  the  crack,  cm 

3 

b  =  co-volume,  cm  /g 

E  =  total  stored  energy,  cal/g 

_  2 
hc  =  local  convective  heat-transfer  coefficient,  cal/cm  -s-K 

F  =  local  convective  heat- transfer  coefficient  over  the  propellant 

p  surface,  cal/cm2-s-K 

h-  =  enthalpy  of  combustion  gas  at  adiabatic  flame  temperature,  cal/g 

J 

k  =  burning  perimeter,  cm 

2 

p  =  static  pressure,  g^/cm 

R  =  specific  gas  constant  for  the  combustion  gases,  gf-cm/g-K 

rb  =  burning  rate  of  the  solid  propellant,  including  the  erosive  burning 

contribution,  cm/s 

T  =  temperature  (without  subscript,  static  gas  temperature),  K 
Tf  =  adiabatic  flame  temperature  of  the  solid  propellant,  K 
Tpt  =  initial  propellant  temperature,  K 

TpS  =  propellant  surface  temperature,  K 

t  =  time,  s 

u  =  gas  velocity,  cm/s 

Vgf  =  velocity  of  propellant  gas  at  the  burning  surface,  cm/s 

x  =  axial  distance  from  propellant  crack  opening,  cm 

xt  =  position  at  the  end  of  crack,  cm 
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NOMENCLATURE  (Cont.) 

y  =  perpendicular  distance  from  the  propellant  surface  into  the  solid,  cm 

o 

a  =  thermal  diffusivity,  cm  /s 

X  ~  thermal  conductivity,  cal/cm-s-K 

p  =  density  (without  subscript,  gas  density),  g/cm 

0  =  angle  measured,  in  a  counterclockwise  direction,  at  the  lower  side 

of  the  propellant,  degree 

~  burning  perimeter,  cm 

Pw  =  wetted  perimeter  of  the  port,  cm 

Subscripts 

i  =  initial  value 

pr  =  solid  propellant  (condensed  phase) 

c  =  rocket  chamber 
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